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Abstract 

We describe recent work of Kim in arXiv:1210.5190 to show that operator 
convex functions associated with quasi-entropies can be used to prove a large 
class of new matrix inequalities in the tri-partite and bi-partite setting by 
taking a judiciously chosen partial trace over all but one of the spaces. We give 
some additional examples in both settings. Furthermore, we observe that the 
equality conditions for all the new inequalities are essentially the same as those 
for strong subadditivity. 
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1 Background 

In a recent paper Kim [IB] showed that for operators on a tensor product space 
Ha ®T~Lb ® 7~Lc one can obtain an interesting new matrix inequality on one space 
He by taking the partial trace Ttab over a quantity for which the full TV abc would 
yield strong subadditivity (SSA) of von Neumann entropy. A key ingredient is the 
inclusion of an additional operator K used in the Wigner-Yanase-Dyson (WYD) 
skew information [30^ 118]. and then choosing a special form for K. As observed at 
the end of [16], Kim's protocol can be applied to a large class of operator convex 
functions, including those associated with the WYD skew information, to produce 
additional matrix inequalities. Although Kim used the recent elegant approach of 
Effros [7], earlier work, going back to Petz [22], based on Araki's relative modular 
operator [3] will suffice. We remark at the end on the different approaches which 
lead to well-known properties used below. 

Kim's result seems remarkable in view of certain well known facts. The concavity 
of the vonNeumann entropy S{p) = —TV p log p is an easy consequence of the much 
stronger operator convexity of g{x) = xlogx. (See, e.g., [BJ and references therein.) 
However, it is also known that the operator h(p, 7) = y^logp — log "n^fp is not 
jointly operator convex, although one has separate operator convexity in the sense 
that the first term is operator convex in p and the second in 7. A trace is needed 
to obtain the joint convexity of relative entropy H (p, 7) = Tr p(log p — log 7) . Thus, 
it is surprising that one can prove an operator inequality whose trace would yield 
something which is an immediate corollary of the joint convexity of relative entropy. 
It should perhaps be emphasized that this does not provide a new proof of SSA; 
rather, the new inequality emerges from the proof of a mild strengthening of SSA 
by the inclusion of an additional operator K as in the WYD setting (from which 
the usual relative entropy can be obtained as a limit) and then making a judicious 
special choice for K in the tripartite setting. 

As he observed, Kim's approach can be used with other convex operator func- 
tions; some examples are worked out below. In addition, we show that in almost all 
cases, the equality conditions are identical to those in [12] for strong subadditivity. 
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2 General theory 



2.1 Some basics 

Let Q denote the class of operator convex functions g(x) on (0,oo) with g(l) = 1. 
For any g 6 Q, one can define a generalized relative entropy [T7], also known as an 
/-divergence [21 H] or "quasi-entropy" [2HI22] as 

H g (K,P,Q) = Tr K*g{L P R£)RQ{K) = Tr K*g(LpR Q r )KQ (1) 

where P, Q > are positive definite matrices and Lp(X) = PX and Rq(X) = XR 
denote left and right multiplication respectively. It is by now well-known that the 
map (P, Q) I—)- H g (K, P, Q) is jointly convex in P, Q for any pair of positive definite 
P, Q and any fixed K. For g(x) = — logx, 

H- logx (K, P, Q) = Tr (k*KQ log Q - KQK* log p) (2) 

which reduces to the usual relative entropy when K = I. 
Whenever g G Q, then g = xg(x~ 1 ) is also in Q and 

Hg(K*,P,Q)=H g (K,Q,P). (3) 
In the case g{x) = — logx above, g = xlogx and 

H xlogx (K,P,Q) = Tr (KK*PlogP- K* PK log Q) (4) 
Because K ^ I is important in what follows, we observe that for any function / 

Tr K*f(L P )KQ = TiK*f{P)KQ 

but 

Tr K*f(R Q )KQ = Tr K*Kf(Q)Q 

which leads to the "sandwiched" expressions KQK* in (|2|) and K*PK in @. This 
is inevitable unless K happens to commute with Q and/or P. 

It is also well-known that joint convexity of H g (P, Q) for fixed K implies mono- 
tonicity under partial traces in the following sense 

H g {I A ® K b c,Pabc,Qabc) > H 9 (Kbc,Pbc,Qbc) (5) 

and that one can not replace I a <8> Kpc by a general Ka.bc ■ The best one can do 
is the minor generalization Va <8> Kpc where Va is unitary [14] . 
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2.2 Main result 



Using the definition (P), we can rewrite ([5]) as 

Tr ABC I a ® K* BC g{Lp ABO RQ 1 ABc )K B c Qabc 

> TT BC K* BC g{L PBC RQ l Bc )K B cQBC (6) 
= Tr abc K* BC g(Lp BC RQ BC )K B cQABC 

where the equality follows from Tic a Qabc = Qbc since there is no explicit depen- 
dence on Ha in the rest of the expression. We also suppress tensor products with 
the identity so that, e.g., it is understood that Pab means Pab ® Ic-, etc. Next, 
consider the special case Pabc — Pab ® Ic to get 

Tvabc I a ® K* BC g{Lp AB R~ l ABc )K BC Qabc 

> Tt AB cKbc9(Lp b Rq bc )KbcQabc (7) 

and choose Kbc = Ib ® Kc- Then K commutes with Pabc — Pab ® Ic so that 
([7|) can be rewritten as 

Tr abc Iab ® K C K C g(Lp AB RQ ABc ) Qabc 

> Ttabc Iab ®K c K c g{Lp B R QBc ) Qabc (8) 

Furthermore, we can choose Kc = |</>c)(</>c| to be a rank one projection so that ([8]) 
becomes 



,Tx A i 



9{Lp ab R q \ bc ) - g^R-^QABC <Pc) > (9) 

where \4>c) is an arbitrary vector in He- This implies that the operator 

TT A B[g{Lp AB R Q \ BC ) - g^R^)] Qabc (10) 
= Ttab 9{L Pab Rq\ bc ) Qabc -Ti B g{Lp B RQ Bc )Q BC > 

is positive semi-definite on He- 

We can choose T-Lp to be one dimensional to obtain bipartite formulas, e.g., 



Tr, 



9( L p ac RQa)- 9(Lpc)] P ac > (11) 



which is also useful. However, we presented the development in the tripartite situ- 
ation because the most important application is to SSA. 
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2.3 Adjoint form 

If one replaces g(x) by g = xg(x _1 ) and interchanges P «->■ Q, one obtains 



^pIRQab) ~ 9(Lpl Rq b )} > (12) 



in which the LHS is formally the adjoint of that in (|1Q|) ; however, since any posi- 
tive semi-definite operator is self-adjoint, (fT2"j) is equivalent to (fTUj) . For additional 
insight into why this is so, recall the basic property that (WX)* = X*W* reverses 
the order and hence, reverses left and right multiplication by self-adjoint operators. 
Thus 

[g{L P R Q l )R Q {X)]* = g(R P L Q l )L Q {X*) 

= L Q R' p l g[{L Q R' p l y l ]Rp{X*) 
= g{L Q R p l )Rp{X*) 

For X = Iab ® \4>c){4 ) c\ = X* the equivalence of (fTUj) and (fT2"j) is then clear. 



3 Specific inequalities 

3.1 Subadditive type 

The choice Qabc — Pabc, Pab = PAB and g{x) = — logx in (fTUj) gives a result 
reminiscent of SSA, i.e., 

Tr^B [ log pabc ~ log Pab - log Pbc + log Pb] PABC > (13a) 

as an operator inequality on He- Using, instead, g(x) = xlogx = xg{x~ l ) and (fT2|) 
with the choices Pabc = Pabc-, Qab = PAB gives the result in the form written by 
Kim in |16| . i.e., 

Tr^B PABC [ log Pabc ~ log PAB + log Pb - log BC ] > (13b) 

which is formally the adjoint of (|13ap . However, as remarked above, these are 
equivalent since a positive semi-definite operator on Tic is necessarily self-adjoint. 
If one uses g = xlogx without the exchange P «-> Q, i.e, with the choice Qabc = 
Pabc, Pab = PAB, one gets 

Tr abPAb[ log pab ~ log pabc ~ log p B + log pbc] > (14) 

in which the simple replacement of pabc by pab on the left appears to reverse the 
usual form of SSA. Note, however, that taking Trc* in (|14p does not yield SSA! 

When Hp is one dimensional, ()13bp becomes 

Tr^ pac [ log Pac ~ log PA ~ log pc] > (15) 
which is an operator version of ordinary subadditivity. 
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3.2 Relative entropy 

The more general choice Qab = 1AB changes (]T3]) to 

Tvab PABC [ log PABC ~ log 1AB ~ log PBC + log Jb] > (16) 

which is naturally associated with the monotonicity of relative entropy under partial 
traces. When T~Lb is one dimensional, this becomes 

TrA Pac [ log Pac ~ log J a ~ log pc*] > (17) 

and for -y A = j^Ia 

- Tr A PAC log Pac + PC log Pc < (log d^) PC (18) 

which gives an upper bound on an operator version of conditional information, 
although this can be negative in the quantum setting. 

Choosing Qabc = Qab <8> Ic m (fl2|) is essential to ensure that it commutes 
with Iab ® Kq- This precludes a proof of the full-fledged operator analogue of 
monotonicity of relative entropy by this method since 

Ti'a pac [ log pac ~ log 7ac - log pc + log 7c] (19) 
is not even Hermitian and, hence, can not be positive semi-definite. 



3.3 WYD inequalities 

The functions g(x) = t ri_ t ) (1 — x t ) and g(x) = t ^i_ t ^ (% — x l ~ l ) generate the WYD 
skew information, for which H g (K, P, P) > and H g (K, P,Q) is jointly operator 
convex in P,Q in the maximal range [—1,2] as observed implicitljo in [3] and ex- 
plicitly by Hasegawa [SI. (See also [H].) Since, as is well known, \img(x) = — logx 

p— >i 

and lim g(x) = xlogx, one can also recover the results of Section [3.11 

Using g(x) in Q with Qabc = PABC, and Pab = IAB gives the inequalities 

t (^-t) TlAB pABC ~ Ttab plj ^c Iab ~ Tt b Pbc - Tr B p]fc 1b > (20) 
for any t G [—1, 2]. Since Tr^B pabc = Tr# psc = PC, this becomes 

[TtabPa'bcIab -Ttbp'bS 1b) > (21) 



t(l-t) 



x Although first Lieb |18j and then Ando [4] obtained the key convexity result for the WYD 
entropy with t G (0, 1), the seemingly innocuous omission of the obvious linear term precludes 
writing their results in the general framework used here. Ando also showed that the concavity of 
Lieb's expression changes to convexity for t £ (1, 2], 
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where it is important to retain the factor mz-x which changes sign at t = 0,1. 
When T~Lb is one-dimensional (|2ip implies 



t(i-t) 



Tr A ^7A > T^DPc' 1 ( 22 ) 



3.4 But Cauchy-Schwarz matrix inequalities are not new 

Using g(x) = (x — l) 2 is equivalent to using x 2 since the linear terms cancel. This 
gives 

T^ab PabQabc ~ r ^ rs PbQbc > (23) 

Since Pab does not depend upon Tic one can use the cyclicity of the trace to rewrite 
this in a more symmetric form as 

Ttab Pab Qabc p ab > Tr B P B Q BC P B (24) 

However, the inequality (|24p is not new; indeed when T~Lb is one-dimensional it 
reduces to something slightly less general than 

Tr A XXcQ-^X A c > X c Qc l X c (25) 

which was provecHin j20] with X—AC arbitrary and Qac positive semi-definite with 
ker Qac Q ker X* AC . Moreover, this is equivalent to the joint operator convexity of 
the map (X, P) i— > X*P~ 1 X also proved in [20] by Lieb and Ruskai, who were 
unaware until 2010 that the latter had been proved much earlier by Kiefer [15] in 
1957. 

However, the slightly modified joint convexity 

(X,P,Q)^TtX* 1 X VtG(0,oo) (26) 

Lp + ttlQ 

does not hold as an operator inequality (even for P = Q) without the trace. Both 
(|25p and (|26p can be proved by very elementary and similar arguments, as shown 
in [20] for the former and for the latter in [25] and the Appendix of [14] . One can 
use (f26j) to prove subadditivity and related inequalities, but the operator inequality 
(I25p does not suffice. This subtle difference makes it even more surprising that 
Kim's method allows one to essentially extract operator inequalities from (|26p in 
the bi-partite ans tri-partite settings. 



In [5D] Lieb and Ruskai proved the slightly more general result that < 
&(X*A~ X) for a completely positive map <E>, of which the partial trace is a special case. 
Later, Choi [5] showed that the hypothesis could be weakened to 2-positivity. The special case 
$(X)*$(X) < was shown earlier for unital maps by Kadison and played an important 

role in Petz's work p2l[2Tl[TT] . 
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3.5 More examples of new inequalities 

Although the functions mentioned above are the most commonly considered, there 
are many more. As shown in [17] , any operator convex function k(x) : (0, oo) i— > 
(0, oo) satisfying the symmetry condition xk(x) = k(x~ 1 ) defines an operator convex 
function g(x) = (1 — x) 2 k{x) with the symmetry property g(x) = x g(x~ l ) = g(x) 
which can be used to define an H g (K, P,Q) as above. The symmetrization 

g{x) + g(x) = — log x + x log x = (x — 1) log x 

yields k(x) = (logx)/(x — 1). However, the symmetrized version yields a less trans- 
parent inequality since one would have the sum of (113aj) and (1141) . 

Several families of functions k have been studied by Petz |23] (who uses f = 1/k 
operator monotone) in the context of monotone Riemmanian metrics and, more 
recently, by Hiai and Kosaki [9j [10] who developed a theory of operator means. A 
fairly comprehensive list is given in [13^ Section 4]. 

The symmetrized version of the primitive example in Section [3.41 is 



\[g{x)+g{x)] = {\ 



x 



1 + x 
2x 



which yields k{x) = (1 + x) /2x. It is well-known |23 [ I17 [ [T3] that the functions k{x) 
satisfy a partial order with 

9 1 -I- r 

< k(x) < — — - (27) 



1 + x ~ K ' ~ 2x 

The smallest element k(x) = 2/(1 + x) is associated with the Bures metric but 
g(LpRg 1 )RQ = (Lp — Rq) 2 (/Lp + Rq) does not seem to yield particularly trans- 
parent inequalities when inserted in (|10l) . 

The function k(x) = x" 1 / 2 also plays a special role in some situations [9[ [TOl 
[T3l 126] . In this case, we can "unsymmetrize" to g(x) = x~ l l 2 — x 1 / 2 and g(x) = 
x(x 1 ^ 2 — x~ 1 / 2 ) to obtain the inequality 

Tr^B 7ab 2 \PABC ~ 1Ab\p 1 Hc ~ Tl B l~p 12 [pBC ~ 1b\Pbc ^ ( 28 ) 
where we used g with P AB = -f AB , Qabc = PABC- 



4 Equality conditions 

It is natural to ask under what conditions equality holds in these inequalities. In 
the case of those related to SSA, i.e., ([l~3"j) and ([I~4"j) . it is easy to see that the 
equality conditions given in [T^ suffice. In the simplest case, Hb = %b' ® 



S 



and pabc = Pab 1 ® Pb"c- The general case is a direct sum of this situation, i.e., 

n B = ® k n k BI ®n k BII and 

PABC = ® P A B> ® P B »C ( 2 9) 

k 

Since a positive semi-definite matrix A > is equal to zero if and only if Tr A = 0, 
it is immediate that (|29p is necessary and sufficient for equality. This also gives 
conditions for equality in (fTQ|) for the other examples with Qabc = Pabc an d 
Pab = Pab and are essentially independent of the function g. This is because 
Nevanlinna's theorem [TJ Section 59, Theorem 2] implies that any operator convex 
function g on (0, oo) with g(l) = has an integral representation of the general 
form 

f°° fix t) 

g(x) = ax + bx 2 + / - v - ; dp g (t) (30) 
Jo x + t 

(The precise representations are written in equivalent, but slightly differently forms 
in several references, including eq. (8.2) in or eq. (17) in |14j or eq. (13) in 
|17] , The details are not relevant here.) Whenever the corresponding measure p g (t) 
is supported on (0, oo), the conditions (j29|) are necessary and sufficient for equality. 
In the approach of |14j the equality conditions arise as the condition for equality in 
the joint convexity in (|26p for all t € (0,oo), which makes the somewhat surprising 
lack of dependence on g transparent. 

For inequalities with more general choices of Pab in (|10j) , as in Sections I3.2L 13.31 
and 13.51 ^ab must also have a similar block representation with p AB , = J^B' f° r 
equality. 

5 Historical remarks 

The proofs of joint convexity in P, Q of functions of the type defined in (pQ) for 
operator convex functions g G Q are based on the relative modular operator Apq = 
LpRg 1 introduced by Araki in a much more general context. The operator L/R 
in Effros's perspective is essentially Apq. Using an integral representation of the 
form ([30]) one can reduce the joint convexity of (P, Q) h-> H g (K, P, Q) to the joint 
convexity of the map in (I26|) which, as mentioned above, can be proved by a very 
elementary argument, as shown in [25] and the Appendix of |14| . 

Uhlmann |28| seems to have been the first to realize that one could take the 
partial trace by integrating unitary conjugations over Haar measure. A pedestrian 
equivalent is to use the discrete Weyl-Heisenberg group (as in, e.g., [Hj), a pro- 
cess sometimes called "twirling" , although other orthogonal unitary bases can also 
be used [27]. Remarkably, Uhlmann [28] also realized that the concavity of the 
map A i— > Tr e K +^°s A could then be used to prove SSA. Without knowing about 
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Uhlmann's work, Lieb found and proved this concavity |18[ Theorem 6 ] which was 
the key ingredient in the two original proofs of SSA presented in [19], both of which 
are different from Uhlmann's. 

Instead of a two-step argument using joint convexity and unitary conjugation, 
one can go directly to the monotonicity under quantum channels, i.e., completely 
positive trace-preserving maps of which the partial trace is a special case. This was 
first done by Petz under the slightly weaker conditional of 2-positivity and a form 
of Jensen's inequality [ITJ ETJ [22]. Another argument was given in [I7j based on 
the integral representation of convex operator functions and an elementary Schwarz 
argument similar to the one in [25] and |14[ Appendix]. 

One advantage of this approach is that one does not need to add the Tr^ to 
the RHS of ([5]) as Kim did in [16] to perform twirling. (See also [14].) One can 
prove ([5]) directly and then use Tr ab ( )ab Qab = Tr abc ( )a Qabc as in the 
last line of 

It is also well-known that one can use a block matrix representation to show that 
monotonicity under partial traces implies convexity, as noted in [18\ fl"9~] [24"] [29] . See 
[Ml ES] for additional background. 
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